Classical and quantum-mechanical treatments of nonsequential double ionization with 

few-cycle laser pulses 
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We address nonsequential double ionization induced by strong, linearly polarized laser fields of 
only a few cycles, considering a physical mechanism in which the second electron is dislodged by 
the inelastic collision of the first electron with its parent ion. The problem is treated classically, 
using an ensemble model, and quantum-mechanically, within the strong-field and uniform saddle- 
point approximations. In the latter case, the results are interpreted in terms of "quantum orbits" , 
which can be related to the trajectories of a classical electron in an electric field. We obtain highly 
asymmetric electron momentum distributions, which strongly depend on the absolute phase, i.e., on 
the phase difference between the pulse envelope and its carrier frequency. Around a particular value 
of this parameter, the distributions shift from the region of positive to that of negative momenta, 
or vice-versa, in a radical fashion. This behavior is investigated in detail for several driving-field 
parameters, and provides a very efficient method for measuring the absolute phase. Both models 
yield very similar distributions, which share the same physical explanation. There exist, however, 
minor discrepancies due to the fact that, beyond the region for which electron-impact ionization 
is classically allowed, the yields from the quantum mechanical computation decay exponentially, 
whereas their classical counterparts vanish. 
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I. INTRODUCTION 

Linearly polarized laser pulses of intensities higher 
than 10 14 W/cm 2 and only a few cycles are of vital impor- 
tance to several areas of physics, being applicable to, for 
instance, solid-state physics 0, high-frequency sources 
, or isolated attosecond pulses . Only the latter ap- 
plication led to a breakthrough in metrology, making it 
possible to trace the motion of bound electrons Q, to 
probe molecular dynamics §, and to control electron 
emission |(J. In this pulse- length regime, the phase dif- 
ference between the pulse envelope and its carrier fre- 
quency, known as "absolute phase" , has a major influ- 
ence on strong-field optical phenomena, such as high- 
order harmonic generation (HHG) Q , or above-threshold 
ionization (ATI) 8]. In particular, the absolute phase af- 
fects, for instance, the harmonic or photoelectron yields, 
the maximal energies in both spectra and the time pro- 
files of ATI and HHG. This is a direct consequence of the 
physical mechanisms governing such phenomena, which 
occur in a subfemtosecond time scale, and for which the 
time dependence of the electric field is important. In 
fact, high-order harmonic generation is the outcome of 
a three-step process in which an electron leaves an atom 
by tunneling ionization at a time t', propagates in the 
continuum, and recombines with its parent ion at a later 
instant t, releasing the energy gained from the field in 
form of high-frequency radiation. A similar mechanism 
is also responsible for ATI, with the main difference that 
the electron either rescatters elastically with its parent 
ion, or reaches the detector without recolliding, originat- 
ing high- or low-energy peaks in the spectra, respectively. 
Such a mechanism has been extensively studied both clas- 
sically [U 0| and quantum-mechanically . 

From the experimental point of view, controlling or 



measuring the absolute phase is a very difficult task [13j • 
This has led to the proposal and experimental realiza- 
tion of schemes for its diagnosis, as, for instance, using 
the asymmetry in ATI photoelectron counts reaching two 
opposite detectors placed in a plane perpendicular to the 
laser field d. 

Another phenomenon whose physical explanation lies 
on a laser-assisted rescattering process is nonsequential 
double ionization (NSDI). In this case, an electron rec- 
ollides inelastically with its parent ion, giving part of its 
kinetic energy to a second electron, which is thus able 
to overcome the second ionization potential and reach 
the continuum. Fingerprints of such a mechanism were 
only revealed very recently, in experiments in which the 
momentum component parallel to the laser field polariza- 
tion could be resolved, either for the doubly charged ion 
[l4|. or for both electrons [l^. Such features, namely a 
doubly-peaked structure in the momentum distributions, 
with maxima at p x \\ — p 2 \\ = ±2^/77p, where Pj\\ (J = 1, 2) 
and Up denote the electron momentum components par- 
allel to the laser field polarization and the ponderomotive 
energy [16|. respectively, are, up to the present date, the 
most striking example of electron-electron correlation in 
the context of atoms in strong laser fields. This fact 
has led not only to further experiments ^tJj but also to 
a considerable theoretical act ivity on the subject, using 
quantum mechanical [Tsl Irfl 12(1, l2lL , semi-classical 
El EHH and classicallHilllainethods. 

Recently, we have shown that NSDI may serve as a 
powerful tool for absolute-phase measurements exploit- 
ing the fact that, for few-cycle driving pulses, inver- 
sion symmetry is broken [29j. Thus, the distributions 
in (j>i||,2?2||) are mainly concentrated in the positive or 
negative momentum regions, changing from one region 
to the other upon a critical phase. Such investigations 
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have been performed classically, considering electrons re- 
leased at times t' uniformly distributed throughout the 
pulse and weighted with the quasi-static tunneling rate 

In this paper, we deal with this problem quantum- 
mechanically, and investigate the existence of a one- 
to-one correspondence with the classical model in Ref. 

Similar studies have been performed for NSDI with 
monochromatic driving fields, with practically identical 
outcomes (2^1 This has shown that, at least in 

the monochromatic case, which is a good approxima- 
tion for the long pulses used in the experiments, intrin- 
sically quantum-mechanical effects such as interference 
processes, or wave-packet spreading, are not important. 
However, it is legitimate to ask the question of whether 
this situation will persist in the few-cycle regime. Indeed, 
it may well be that interference and wave-packet spread- 
ing play a more important role in this latter case. Addi- 
tionally, it is not clear whether the quasi-static tunneling 
rate considered in the classical model remains valid for 
few-cycle driving pulses. In fact, this has been recently 
put into question, with the derivation of a non-adiabatic 
rate [3l|. Finally, it is worth checking whether asym- 
metric distributions and the critical phase also occur in 
a quantum-mechanical context, and, in case they do, to 
understand the physics behind. 

In particular, we address the above-stated questions 
using a S-Matrix formalism, within the Strong-Field Ap- 
proximation (SFA) [33, We consider the simplest 
type of rescattering, namely electron-impact ionization, 
and treat the problem in terms of the so-called "quantum 
orbits" [34|, which appear in the context of saddle-point 
approximations. Specifically, we use a uniform approxi- 
mation whose only validity requirement is that the orbits 
in question occur in pairs, which is in general the case for 
laser-assisted rescattering phenomena. This method has 
been previously applied to NSDI in monochromatic driv- 
ing fields, in order to analyze the influence of the types 
of interaction and final-state electron correlation on the 
yields [2{J [2(| . Apart from considerably simplifying the 
computations involved, as compared to other theoretical 
methods 0, 0, |2(j , the quantum-orbit approach pro- 
vides additional physical insight, in terms of a space-time 
picture. In fact, the quantum orbits are closely related 
to the orbits of a classical electron in an external laser 
field. Hence, in several situations, it is possible to draw 
a parallel betweeen our quantum-mechanical treatment 
and the previous classical considerations [29| . discussing 
their similarities and differences. In the following, we 
study the physical mechanisms responsible for the critical 
phase within a quantum-mechanical framework, concen- 
trating on the main differences from the classical picture 
and from the monochromatic-driving field case. 

The paper is organized as follows: In the next sec- 
tion (Sec. FUjl. we provide the necessary theoretical 
background, presenting the transition amplitude in the 
strong-field and uniform approximations. Subsequently, 
in Sec. IIIII we present differential electron momentum 



distributions for various absolute phases, discussing the 
main features obtained in terms of quantum orbits. The 
quantum mechanical results are then compared to a clas- 
sical ensemble computation which is either the same as 
in 29], or slightly modified with respect to it (Sec. HVfl . 
Finally, in Sec. ]V] we summarize the paper and state our 
conclusions. 



II. BACKGROUND 

The transition amplitude of the laser-assisted inelastic 
rescattering process responsible for NSDI, in the Strong- 
Field Approximation 32, 33] , is given by 
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where V, U^\t,t'), U^ v> {t,t'), and V 12 denote the 
atomic binding potential, the field-free and the Volkov 
time evolution operators acting on the n-th (n = 1,2) 
electron, and the interaction through which the sec- 
ond electron is freed by the first, respectively. Eq. 
(Q) expresses the following physical process: Initially, 
both electrons are bound, and the atom is in the 
ground state, which is approximated by \ipo(t')) = 

ip^\t')\(i.e., product state of one-electron 
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ground states), with ^o"^*')) = e 4 ^ ™ 1 *' 

time t' , the first electron is released through tunneling 
ionization, whereas the second electron remains bound. 
Subsequently, the first electron propagates in the contin- 
uum from t' to t, gaining energy from the field. At this 
latter time, it collides inelastically with its parent ion, 
dislodging the second electron. The final electron state is 
then chosen as the product state of one-electron Volkov 



waves, 



4l ,P2 {t) 



4V 



(£)), where pi,p 2 



are the final electron momenta (for studies of correlated 
two-electron final states see, e.g., Refs. ^8j, [2^ and 
[2^]). In Eq. Q), the interaction with the ionic potential 
is not taken into account. In our computations, we use 
the length gauge and atomic units. 

Expanding U^ v \t,t') in terms of Volkov states, Eq. 
reads 

dt / dt' / d 3 kV pk Vwexp[iS(t,t',pi,-k)], 

-OO J — OO J 
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with the action 

S(t,t', Pl ,k) 



(2) 



i=i 



[P, + A(r)] 2 dr 



i^[k + A(r)] 2 dr 

|Bbi|t' + |^|t, 



(3) 
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where A(t), Pj(j — l,2),k and \Eo n \ denote the vector 
potential, the final momenta of both electrons, the in- 
termediate momentum of the first electron and the ion- 
ization potentials, respectively. All the influence of the 
binding potential V and of the electron-electron interac- 
tion V12 is included in the form factors 

Vp k =<p 2 + A(t), Pl +A(i)|yi 2 |k + A(i),4 2) > (4) 
and 

y k0 -<k + A(t')l^|^ 1) > • (5) 

In this paper, we consider a contact-type interaction 

^12=5(^-^)^2), (6) 

which yields very good agreement with experimental data 
within the context of NSDI in monochromatic driving 
fields [23|,|25j. In this case, the form factors Vpk, Vko are 
constant and the SFA transition amplitude can be solved 
analytically up to one quadrature. For other types of 
potentials, this is only possible by evaluating multiple 
integrals numerically. 

For low enough frequencies and high enough laser in- 
tensities, Eq. J2| can be solved to a good approximation 
by the steepest descent method. Thus, we must deter- 
mine k, t' and t, such that S(t, t' , p^, k) is stationary, i.e., 
its partial derivatives with respect to these parameters 
vanish. This condition yields 

[k + A(t')} 2 = -2\E 01 \, (7) 



2 

5>,+A(t)]Mk + A(t)] 2 -2|i; 02 | (8) 
i=i 

and 

J dr [k + A(r)] - 0. (9) 

Eq. and (jSJ) give the energy conservation at the start 
and rescattering times, respectively, while Eq. con- 
straints the intermediate momentum of the first electron 
so that it returns to its parent ion. For vanishing \Eqi\, 
the classical equations of motion of both electrons in the 
external field are obtained. For non-zero |2?oi|j Eq. (JTJ) 
expresses tunneling ionization at t' , and has no real so- 
lution. Physically, this means that this process is not 
classically allowed. This results in complex variables t', t 
and k, which always occur in pairs. The real parts of such 
variables are directly related to a longer and a shorter or- 
bit of a classical electron in an electric field. The longer 
orbit can be associated to the so-called "slow-down colli- 
sions", which have recently been discussed in the litera- 
ture |22ll24. 28] . The imaginary parts determine to which 
extent electron-impact ionization is allowed or forbidden, 
both within and beyond the boundaries of the classically 



allowed energy region. In this latter domain, one of the 
orbits leads to exponentially decaying contributions in 
the transition amplitude @, which cause cutoffs in the 
distributions, while the remaining orbit starts to yield 
diverging contributions, and must be discarded. 

Eq. (jSJ) can also be written in terms of the momentum 
components parallel and perpendicular to the laser-field 
polarization, denoted by pj\\ and Pj±(j = 1,2), respec- 
tively. In this case, for constant transverse momenta, one 
obtains the equation 

2 2 

f^PiU + m? = [k + A(t)] 2 -2\E 02 \ -J2p 2 3± , (10) 

describing a circle in the p^u, p 2 \\ plane, whose radius de- 
pends on the kinetic energy E lct (t) — [k+ A(i)] 2 /2 of the 
first electron upon return, and on the effective binding 

energy |£ 02 | = \E 02 \ + £ p 2 ± /2. If E Tet (t) < \E 02 \, this 

radius collapses and electron-impact ionization becomes 
classically forbidden. This means that there is not only 
a maximal, but also a minimal classically allowed energy, 
and the resulting yields exhibit two cutoffs, or no cutoff 
at all. This is a major difference with respect to high- 
order harmonic generation or above-threshold ionization, 
for which only maximal classically allowed energies exist. 

In the standard saddle point method, the action © 
is expanded quadratically around the saddle points, and 
the transition amplitude @ is approximated by 

m( spa) = J2A s exp(iS s ), (11) 

S 

S s = S v (t s ,t' s ,K), (12) 

As = (2tu) 5 / 2 . Vpk * Fk *° , (13) 
^detS£(M',k)| s 

where the index s runs over the relevant saddle points, 
and Sp denotes the five-dimensional matrix of the second 
derivatives of the action with respect to t, t' and k. In 
practice, we first determine k(i,i') as a function of the 
other variables, inserting this in the action, and take 

A s = (2m) 5 / 2 V p ^o (14) 

(i'-t)3/2^detS£(i,tf)| s 

so that the computation of the determinant is simplified. 

The above-stated saddle-point approximation is only 
applicable for well-isolated saddle-points. This does not 
hold near the boundaries of the classically allowed region, 
where the pairs of saddles nearly coalesce. Furthermore, 
beyond such boundaries, one of the saddles yields diverg- 
ing results, and must be discarded. This leads to cusps 
in the yield which are particularly problematic for non- 
sequential double ionization. A detailed analysis of this 
problem is given in [24|. 

Such artifacts can be eliminated by using a more gen- 
eral, uniform approximation [37l | , whose only validity re- 
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quirement is that the saddles occur in pairs. This ap- 
proximation has been successfully applied in the context 
of above-threshold ionization j^fj and nonsequential dou- 
ble ionization jHHl 

Within this improved approximation, in the classically 
allowed region, the transition amplitude for a pair of tra- 
jectories i and j is given by 



M, 



i+j = y / 2nAS/3 exp(iS + in/ A) 

x {A[J 1/3 (AS) + J_ 1/3 (A5)] 
+ AA[J 2/3 (AS) - J_ 2/3 (A5)]}, 
AS = {Si - Sj)/2, S = ($ + Sj)/2, 
AA = (Ai-iAj)/2, A=(iA i -A j )/2. 



(15) 



The saddle-point approximation is recovered for large 
values of AS, using the asymptotic behavior 



J±u{z) ~ — 

7TZ 



1/2 



cos(z =F - tt/4) (16) 



of the Bessel functions for large z. One should note that 
the uniform approximation considers the collective con- 
tribution of a pair of saddle points, instead of, as in the 
former case, taking them into account individually. 

In the classically forbidden region, one of the saddles 
must be discarded. For that purpose, a branch of the 
Bessel function must be chosen in such a way that the 
approximation exhibits a smooth functional behavior at 
the Stokes transition [3{| , given by 

Re5p(t < ,4k i ) = ReSp(t i ,tJ.,Itf). (17) 
Beyond the Stokes transition, 



M i+j = y/2iAS/-Kexp(iS) 

x [AK l/3 (-iAS) + iAAK 2/3 (-iAS)] . (18) 

The saddle-point approximation is, again, recovered us- 
ing the asymptotic expansion 



K, 



/ 7T \ V 2 

.(*) ~ (-] cxp(-z) 



(19) 



for large z. Inserting Eq. I|19|) into (|18|l . it is easy to show 
that only one saddle contributes to the saddle-point ap- 
proximation in this energy region. Equations (|15J) and 
(|18|l should be matched at the Stokes transitions, whose 
energy positions roughly coincide with the boundary be- 
tween the classically allowed and forbidden energy re- 
gions. 

In the following, we take the few-cycle pulse E(f) = 
-dA(t)/dt, with 

A(i) = A Q exp[-4(wi - irn) 2 / '(™) 2 ] sin [a;* + <f>]e x , (20) 

where n, u>, A and cj> denote its number of cycles, fre- 
quency, amplitude and absolute phase, respectively. We 
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FIG. 1: Differential electron momentum distributions com- 
puted for neon (|-Eoi| = 0.79 a.u. and \Eoa\ = 1.51 a.u.) sub- 
ject to a four-cycle pulse (n = 4) of freqency uj = 0.057 a.u. 
and various intensities and absolute phases. The upper, mid- 
dle and lower panels correspond to I = 4 x 10 14 W/cm 2 ([/ p = 
0.879 a.u), I = 5.5 x 10 14 W/cm 2 ([/ p = 1.2 a.u), and I = 
8 x 10 14 W/cm 2 (l/ p = 1.758 a.u), respectively. The absolute 
phases are given as follows: Panels (a), (e) and (i): <j> = 0.87r; 
panels (b), (f) and (j): <j> — 0.97r; panels (c), (g) and (k): 
cj> — Itt; and panels (d), (h) and (1): <f> = 1.1-7T. 



then find the start and return times such that the saddle- 
point equations are fulfilled, and use such times to com- 
pute the yields, which are given by 



r(pi||,P2||) = J d 2 p lx J d 2 p 2 x\M\ 



(21) 



where M is given by Eq. J2J) within the uniform approx- 
imation. 



III. QUANTUM-ORBIT ANALYSIS 

In Fig. 1 we present the momentum distributions com- 
puted using the above-discussed method, for various ab- 
solute phases, in form of contour plots in the (f>i||,P2||) 
plane. We choose the atomic species to be neon, for 
which electron- imp act ionization is the dominant physi- 
cal mechanism [38| . In general, such distributions exhibit 
circular shapes, centered at particular momenta along 
Pi || = f 2 1| = P||) and, in contrast to the monochromatic- 
field case, are no longer symmetric with respect to 
(Pi\\,P2\\) (-Pi||,-P2||)< This symmetry breaking is ex- 
pected, since the relation A(t) = —A(t ± T/2), and thus 
\M(t,t',p lh p n )\ = \M(t ± T/2,t' ± T/2, -p 2 ||)|, 
which was true for monochromatic driving fields, does not 
hold. The circular shapes are typical for the contact-type 
interaction, and are also observed in the monochromatic 
case. 

Depending on the phase, the yields are mainly concen- 
trated either in the regions of positive or negative parallel 
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FIG. 2: Real and imaginary parts of the start and return times 
for the orbits 1 and 2 as functions of the parallel momentum jjm 
along the diagonal pi\\ = Pay , computed for a four-cycle pulse 
(n = 4) of absolute phase <f> = 0.5n. The atomic parameters 
were taken as |-Eoi — 0.79 a.u. and |Soa| = 1.51 a.u. and 
correspond to neon, while the field intensity and frequency 
were chosen as I = 5.5 x 10 14 W/cm 2 (i7j, = 1.2 a.u) and ui = 
0.057 a.u, respectively. The numbers in the figure denote the 
transverse momenta (pi_i_,paj_) in units of *JU P . The shorter 
and longer orbits in each pair correspond to the solid and 
dashed lines, respectively. 
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FIG. 3: The same as in the previous figures for the orbits 3 
and 4. 



momenta. For instance, in the figure, initially, the par- 
allel momenta of both electrons are essentially positive 
(Figs. l.(a), l.(e) and l.(i)). As the phase increases, 
contributions from negative momenta are also present, 
becoming more and more significant, until the distribu- 
tions are almost entirely shifted from the positive to the 
negative momentum region (c.f. Figs. l.(d), l.(h) and 
1.(1)). This process occurs for different intervals of abso- 
lute phases, depending on the peak intensity of the driv- 
ing field. For the specific example presented, the higher 
the intensity is, the earlier the momenta start to change 
sign. This phase dependence is very similar to that in 



[29j , obtained within a classical framework. Thereby, this 
behavior was traced back to sets of electron trajectories, 
whose relevance was determined by the phase space, and 
by the rate with which the first electron was ejected in 
the continuum. A critical phase was related to a change 
in the dominant pair of orbits, which had a huge reper- 
cussion in the distributions. This phase was also shifted 
towards smaller absolute values with increasing driving- 
field intensity. 

Subsequently, we analyze both the asymmetry and the 
critical phase in terms of pairs of quantum orbits, which 
are classified as according to increasing start 

times and absolute phases. We consider only relatively 
short orbits so that t — t' <T, where T = 2tt/uj denotes 
the field cycle. Longer orbits yield negligible contribu- 
tions due to wave-packet spreading. 

Fig. 2 shows one of such pairs for t' near T, which 



we specifically denote (1, 2) 



(0.5tt) 



We consider the inter- 



mediate intensity on Fig. 1, for which the phase cho- 
sen yields positive parallel momenta, along the diagonal 
Pi\\ — P2\\ — P\\- Panels (a) and (b) display Re[wt'] and 
Re[u;i] as functions of py, which can be associated to the 
times obtained by solving the classical equations of mo- 
tion of two electrons colliding inelastically in a laser field. 
There exists always a longer and a shorter orbit for the 
electron, which nearly coalesce near two particular mo- 
menta. Such momenta correspond to Stokes transitions 
[Eq. I|17l) ]. which, for high enough intensities, roughly co- 
incide with the minimal and maximal classically allowed 
momenta [39| . These two specific momenta delimit a 
region, which is most extense for = (j = 1,2). 
As the transverse momenta increase, the effective second 
ionization potential |2?o2| also becomes larger until this 
region collapses. An interesting feature is that, between 
the Stokes transitions, the real parts of the rescattering 
and start times are centered around a particular value 
of pi: , which correspond to the peak of the momentum 
distributions. For few-cycle pulses, this center depends 
on the pair of orbits, as well as on the absolute phase. 
For monochromatic driving fields, it lies at pu = ±2^/[7p 
0. 

The remaining panels depict the imaginary parts of 
such times, which provide in some sense a measure for 
a process being classically allowed or forbidden. Indeed, 
they determine whether the transition amplitudes J5J in- 
crease or decrease exponentially, or how relevant the con- 
tributions from particular sets of orbits are. These imag- 
inary parts noticeably increase at and beyond the Stokes 
transitions, and remain practically constant in the mo- 
mentum region inbetween. This suggests that, in this 
region, electron-impact ionization is either classically al- 
lowed, or at least much more probable to occur. Ad- 
ditionaly, whereas Im[cji] almost vanishes in this region, 
Imjwt'] has a nearly constant and nonvanishing value (c.f. 
panels (c)). This is due to the fact that tunneling per se 
is classically forbidden. Indeed, the larger this value is, 
the smaller is the probability that this process takes place 
at all. 
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FIG. 4: The same as in the previous figures for the orbits 5 
and 6. 



Two additional pairs of orbits, for which 1.5T < t 1 < 
2T, are displayed in Figs. 3 and 4. At these times, the 
pulse l|2U|) is closer to its peak intensity. In such figures, 
there exist extense regions between the Stokes transi- 
tions, in which Re[wi'] and Re[wi] practically coincide 
with the start and return times obtained within a clas- 
sical framework, and in which |Im[u;£]| are vanishingly 
small. Such features are a clear evidence that, in this 
case, electron-impact ionization is classically allowed. 

Another noteworthy feature is that, in the classically 
allowed region, |Im[w£']| has much smaller values than 
those in Fig. 2. Physically, this means that the first 
electron left with a larger tunneling probability at t' , in 
comparison to the orbits (1.2)(o.57r). Furthermore, the 
fact that this region is extense shows that that the ki- 
netic energy of the first electron upon return is larger for 



(3,4) 



(0.5tt) 



and (5,6)( .5 7r ) than for (1.2) 



(0.5tt)- 



For that 



reason, the contributions from (3,4)(o.5,r) and (5,6)(o.5 W ) 
to the total yield should be more relevant than those from 

( 1 - 2 )(0.5ti-)- 

This is confirmed by Fig. 5, which depicts the yields 
computed from each pair of the above-discussed orbits, 
along the diagonal p^n = p 2 || = P\\- I 11 Fig. 5(a), the 
contributions from (5,6)(o.5 ff ) are at least two orders of 
magnitude larger than those from the remaining pairs, 
so that the distributions will be concentrated in the first 
quadrant of the (f>i||,P2||) plane. Hence, for practical 
purposes, the remaining contributions can be neglected. 
They are, however, very useful for the physical under- 
standing of the problem. 

The second most prominent contributions come from 
(3,4)( 5 ir). This is expected, since, for these orbits, there 
is a relatively large probablity that the first electron tun- 
nels out, as well as a large momentum region for which 
electron-impact ionization is allowed. 

Additional contributions come from the orbits 
(1, 2)(o.5tt) and (7, 8)(o.5 ff ). The latter set of orbits is not 
displayed in the previous figures, due to the fact that, 
in this case, there are no Stokes transitions, i.e., elec- 
tron impact ionization is forbidden throughout. Interest- 




FIG. 5: Individual contributions of the four most relevant 
pairs of orbits to the NSDI yield, for p^i = p 2 u = py, for pon- 
deromotive energy U p = 1.2 a.u and absolute phases <j> = 0.5n, 
cj> — 0.8-7T and <j> = l.lir (panels (a), (b) and (c), respectively). 
The remaining parameters are the same as in the previous 
figures. The curves have been normalized to the maximum 
of the most relevant contributions. Specifically, in panels (b) 
and (c), the contributions from (1,2) are smaller than the 
range of orders of magnitude displayed. 



ingly, the contributions from all pairs of orbits discussed 
above, including (7, 8)(o.5 ff ), are several orders of magni- 
tude larger than those from the pair (1, 2)(o.5,r). This is 
due to the fact that, for (7, 8)(o.5 ff ) the tunneling prob- 
ability for the first electron is considerably larger than 
for (1, 2)( 57r ). From the technical viewpoint, it is worth 
mentioning that, for (7, 8)(o.5 V ), the yield has been com- 
puted by using Eq. I|ll[) . and taking the orbit for which 
this expression is exponentially decaying. 

For other absolute phases, there may be other sets of 
orbits whose contributions may compete with or even 
overwhelm those from (5,6). This is in fact the case in 
Figs. 5(b), and 5(c), for <\> = 0.87T and <j> — 1.1-k, respec- 
tively. Such phases, as well as the remaining parameters, 
are the same as in Figs. 1(d) and 1(f), corresponding to 
the beginning and to the end of a shift in the momentum 
distributions. 

In Fig. 5(b), one clearly sees that the second most rele- 
vant pair of orbits is no longer (3, 4)(o.8tt), but (7, 8)(o.8ir)- 
The contributions from such orbits now are only one or- 
der of magnitude smaller than those from (5, 6)( .8ir)- 
Consequently, there are also small, but not negligible 
contributions in negative momentum regions. This is in 
agreement with Fig. 1(d), for which there is a small spot 
in the third quadrant of the (pi||,P2||) plane, in addition 
to the dominant contributions in the first quadrant. For 
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FIG. 7: Time-dependent electric field, for various absolute 
phases. The peak-field times are marked with arrows and the 
corresponding pairs of orbits are indicated by the numbers in 
the figure. The remaining parameters are the same as in the 
previous figure. 



FIG. 6: Real (panel (a)) and imaginary (panel (b)) parts 
of the start times for the orbits (3,4), (5,6) and (7,8) along 
Pill — P211 an d for vanishing transverse momenta, indicated 
by dotted, dashed and solid lines, respectively. We consider 
neon subject to a four-cycle (n = 4) pulse of intensity / = 
5.5 x 10 14 W/cm 2 , various absolute phases and frequency cu = 
0.057 a.u. Specifically for <j> — 0.5-/T, electron-impact ionization 
is classically forbidden for the orbits (7,8). 

larger phases, as for instance <j> = l.ln (Fig. 5(c)), the 
contributions from (7, 8) become even more relevant than 
those from (5,6). Hence, the distributions are shifted 
from the first to the third quadrant, in accordance to 
Fig. 1(h). 

In Fig. 6, we systematically analyze the dependence of 
the real and imaginary parts of the times t' on the abso- 
lute phase, for the three most relevant sets of orbits. As 
an overall feature, for the pair (3,4), |Im[a;f']| is larger 
than for the remaining two pairs. Physically, this means 
that tunneling is less probable in this case. In general, 
this makes the contributions from this pair to the yield 
less relevant than those from the other two sets. An ex- 
ception, however, occurs for <f> = 0.5n, due to the fact 
that, in this case, electron-impact ionization is forbidden 
for (7, 8). This leads to exponentially decaying contribu- 
tions for this pair, which are overwhelmed by those from 
(3,4). Still, even in this particular case, these latter are 
two order of magnitude smaller than the yield from the 
orbits (5,6) (c.f. Fig. 5(a)). 

For high driving-field intensities, as those in the lower 
panels in Fig. 1, there may exist minor, though not neg- 
ligible, contributions from (3,4). In general, however, 
these contributions are vanishingly small. An interesting 
feature is that, with increasing absolute phases, this set of 
orbits gradually loses relevance, since |Im[u;t']| increases. 
This can be explicitly seen in Fig. 5, where the yield from 



(3,4) decreases in at least three orders of magnitude as 
4> increases. 

The other two pairs of orbits, (5,6) and (7,8) are, in 
fact, far more important to the yield. Thereby, three dis- 
tinct behaviors can be identified. Below <f> = 0.8ir, the 
distributions are mainly determined by the orbits (5,6), 
whose contributions lie in the region of pn > 0. Around 
this phase, the pair (7, 8)( .8ir) comes into play. Indeed, 
although this pair does not delimit a large classically al- 
lowed region, the imaginary parts of the corresponding 
start times are comparable to or smaller than those of 
(5,6). Thus, the contributions from both pairs start to 
compete, and the distributions spread over both positive 
and negative momentum regions. As the phase increases, 
|Im[wt']| gradually decreases for this latter set, until the 
negative parallel momenta are favored. This explains the 
features in Fig. 1. One may refer to <f> = 0.8tt as a critical 
phase 4> c , since it marks a change in the sets of dominant 
orbits. 

The pulse profile (Fig. 7), together with the real parts 
of the tunneling times, allows an intuitive interpretation 
of the above statements. Below <p c = 0.8n, the peak 
value of the pulse is near Re[wt'] for the orbits (5,6), 
whose contributions then dominate. Thus, the distribu- 
tions essentially concentrate in the positive momentum 
region. Around this phase, this picture starts to change, 
and there are two sets of orbits, namely (5, 6) and (7, 8), 
for which the instantaneous electric fields at the tunnel- 
ing times are comparable. This situation persists within 
a phase interval until, finally, the absolute maximum of 
the field corresponds to the latter set of orbits, so that the 
momenta are mainly negative. An interesting situation 
occurs for (f> = 0.57T, for which, in principle, there exists 
two sets of times near which the electric field exhibits 
comparable maxima, corresponding to (3,4) and (7,8). 
For low intensities, rescattering is not allowed for the or- 
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FIG. 8: Real (panel (a)) and imaginary (panel (b)) parts of 
the start times for the orbits (5,6), (7,8), indicated by dashed 
and solid lines, respectively, for absolute phase </> = 0.8-7T and 
various laser intensities, along pxy = p 2 | and vanishing trans- 
verse momenta. The remaining parameters are the same as 
in the previous figure. 



bits (7,8), and the contributions to the yield are absent. 
For high intensities, however, electron-impact ionization 
is already allowed in this case, so that the critical phase 
is shifted towards smaller values. 

In Fig. 8, we explicitly show how electron-impact ion- 
ization becomes classically allowed or or forbidden upon 
a change in the driving-field intensity, for <f> c = 0.8-7T. In 
fact, as the intensity is decreased, this process ceases to 
be allowed for the orbits (7,8), whose contributions lie 
near the peak of the pulse (c.f. Fig. 7). As a direct 
consequence, the critical phase is shifted towards larger 
absolute values as the intensity decreases, as shown in 
Fig. 1. 



IV. COMPARISON WITH CLASSICAL MODELS 

In this section, we will recall and apply the classical 
model used in 29] , performing a direct comparison with 
the results of the S'-Matrix computation. We consider 
an electron ensemble subject to the few-cycle pulse (|2Up. 
which are released from the origin of the coordinate sys- 
tem with initial vanishing drift velocities, i.e., 



p + A(t') = 0. 



(22) 



The varying parameters are the tunneling times i', uni- 
formly distributed throughout the pulse, and the quasi- 



R(t>) 



exp 



-2(2| J B 01 |) 3 / 2 /3| J B(^)ll /\E(t% (23) 



with which the electron counts are weighted. Some of 
these electrons subsequently return to the origin and 
free a second electron ensemble at later times t through 
electron-impact ionization. Their return and rescattering 
conditions are given by Eqs. 10 and JSJ, respectively. 
The yields are then given by 



T ~ J dt'R(t')S ^Er, 



*(*)-E 



(Pi + A(t)) s 



I -Bi 
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where E let (t) is the kinetic energy of the electron upon 
return and the argument in the S function gives the en- 
ergy conservation at t. This model is discussed in more 
detail in H^. 

Fig. 9 presents the outcome of the classical compu- 
tation, for the same parameters as in Fig. 1. Inter- 
estingly, both figures are very similar. Indeed, there 
exist only minor differences near the boundaries of the 
classically allowed region, which occur for high inten- 
sities, as displayed in the lower panels of Figs. 1 and 
9. Such differences are due to the fact that the yield 
from the quantum-mechanical computation is exponen- 
tially decaying in the region for which electron-impact 
ionization is classically forbidden, whereas the distribu- 
tion 124|) immediately vanishes. Such discrepancies were 
also present in the monochromatic case. 

Further discrepancies occur for specific phases, in pan- 
els (c) and (e) of both figures, and show that the sign 
reversal starts to take place for slightly smaller abso- 
lute phases in the quantum-mechanical case. This effect 
can be understood if one keeps in mind that the criti- 
cal phase indicates a change in the dominant set of tra- 
jectories, and that pre-requisites for this dominance are 
a large tunneling probability for the first electron and 
electron-impact ionization being classically allowed. In 
the particular example provided in the panels (c) and (e) 
of Figs. 1 and 9, this process has just become allowed 
for the orbits (7, 8), within a small momentum region. In 
the quantum- mechanical case (Fig. 1), there exist con- 
tributions to the yield from near the boundary of such a 
region, even if rescattering is forbidden, whereas in the 
classical model (Fig. 9) such contributions vanish. Ob- 
viously, such discrepancies are absent in the symmetric 
momentum distributions obtained in the monochromatic 
case. 

For monochromatic driving fields, a very good agree- 
ment has also been reported in previous publications 
[2^,|2^|. It is not obvious, however, that this would re- 
main true in the few-cycle regime. In fact, apart from 
yielding distributions which immediately vanish at the 
boundaries of the classically allowed region, the classical 
model does not take into account several effects which 
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FIG. 9: Differential electron momentum distributions for the 
same parameters as in Fig.l, computed with the classical 
model discussed in Sec. IV. 

are present in the quantum-mechanical computation. Ex- 
amples of such effects are the quantum interference be- 
tween different possible paths for the returning electron, 
or the spread of the electron wave packet. Furthermore, 
the classical model considers an additional approxima- 
tion with respect to the ^-Matrix computation, namely 
Eq. J2BJ), which is a quasi-static, cycle- averaged tunnel- 
ing rate. This rate is a key ingredient in the distributions 
given by Eq. (|24[l . and vital for the phase-dependence 
observed in Fig. 9 (c.f. Fig. 10 and discussions in [2^|). 
However, its validity may be limited or even questionable 
in the few-cycle regime. Indeed, recently, a non-adiabatic 
rate that should be more accurate in this case has been 
derived [3]| . In general, the non-adiabatic rate tends to 
broaden the time range for which the relevance of sets 
of orbits persists. We verified, however, that this rate 
does not modify the yields for the parameter range in 
question. Discrepancies between the quasi-static and the 
non-adiabatic rate occur only when the Keldysh param- 
eter 7 = y/\Eoi\/(2Up) is much larger than unity. In this 
case, the driving-field intensity would be far too low for 
the rescattering process discussed in this paper to be al- 
lowed, and thus for the classical model to be applicable. 
For measurements of NSDI yields below the threshold 
see, e.g., Ref. [loj. 

The subsequent picture (Fig. 10) is the classical coun- 
terpart of Fig. 6, providing an interpretation of Fig. 
9 in terms of the interplay between the ionization rate 
()23J) and the phase space. Thereby, as in the quantum- 
mechanical case, we consider parallel momenta along the 
diagonal p^ = p 2 \\ = p\\ and vanishing transverse mo- 
menta. The main contributions to the yield come from 
pairs of trajectories for which the tunneling probability 
is large, and for which electron-impact ionization is clas- 
sically allowed. In the figure, according to such criteria, 
it is possible to identify two sets of relevant orbits, cor- 
responding to electrons ejected at 9 < u>t' < 12, with 
positive parallel momenta, and to electrons released at 



FIG. 10: Quasi-static tunneling rate t[2'A} fnanel (a)), together 
with the classically allowed momenta computed with the clas- 
sical model (panel (b)), as functions of the tunneling times, 
for the same parameters as in Fig. 6. 

12 < u>t' < 15, with negative parallel momenta. 

According to the absolute phase, similarly to the 
quantum-mechanical case, there exist three distinct types 
of behavior. For <j> < 0.8ir, apart from the fact that 
the classically allowed region is almost vanishing for the 
orbits (7,8), the quasi-static rate (|23[) is considerably 
larger for the former set of electrons, so that the distribu- 
tions are concentrated in the positive momentum regions. 
Around <f> = 0.87T, such tunneling rates are comparable 
for both sets of orbits, resulting in nonvanishing yields for 
positive and negative parallel momenta. Finally, as the 
phase increases, tunneling is more prominent for the lat- 
ter set of orbits, and the distributions gradually change 
towards negative momenta. As in the previous section, 
the critical phase marks a change in the dominant orbits. 
This is in agreement with Fig. 6 and with the results in 

The ratios between the individual contributions of such 
pairs of orbits, as well as the momenta for which their 
maxima occur, displayed in Fig. 11, are also in very 
good agreement with its quantum-mechanical counter- 
part (Fig. 5). However, the distributions computed with 
the classical model are narrower than those obtained us- 
ing the quantum-mechanical computation. This is once 
more related to the fact that, in a classical framework, 
the contributions from the forbidden momentum region 
can not be taken into account. 

Another interesting feature, which is displayed in Fig. 
12, is the existence of other critical phases. For instance, 
around </> = 0.17T, there is a transition in the momentum 
distributions from the third to the first quadrant in the 
(Pi||iP2||) plane, i.e., exactly in the opposite direction as 
the transition in Figs. 1 and 9. Furthermore, the upper 
and lower panels in Fig. 12 look exactly the same, if the 
first and the third quadrant are interchanged. This shows 
that there is a symmetry in the momentum distributions, 
which is due to the fact that A(t, <fi) = — A(t, <fi ± tt), so 
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FIG. 11: Individual contributions to the momentum distribu- 
tions for the same pairs of orbits and for the same parameters 
as in Fig. 5, computed with the classical model. Parts (a), 
(b) and (c) correspond to <j> = O.bn, <f> = 0.87T and (f> = l.lir, 
respectively. The curves have been normalized to the peak 
value of the largest curve. The contributions from the orbits 
for which electron-impact ionization is classically forbidden 
are absent. 
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FIG. 12: Differential electron momentum distributions com- 
puted with the classical model in Sec. IV for neon sub- 
ject to a four-cycle pulse (n = 4) of approximate intensity 
7 = 5.5 X 10 14 W/cm 2 (Up=1.2), frequency u> = 0.057 a.u. 
and absolute phases <j> = O.ln, <j> = 0.57T, and (j> — 0.9n (pan- 
els (a) to (c), respectively), and <f> = 1.17T, <f> = 1.5-7T, and 
cj> — 1.97T (panels (d) to (f), respectively). 



that \M( Pl[l ,p 2ll ,4>)\ = \M(-p lh -p n ,(P±Tr)\. 



V. CONCLUSIONS 

The studies performed in this paper clearly show that 
nonsequential double ionization with few-cycle pulses is 
a powerful tool for absolute-phase measurements. More 
specifically, the yields, as functions of the electron mo- 
mentum components (pi||,P2||) parallel to the laser-field 
polarization, are mainly concentrated either in the posi- 
tive or negative momentum region, depending on the ab- 
solute phase in question. Around a critical phase, such 
distributions start to shift from one momentum region 
to the other, until, as the phase increases, complete sign 
reversal in the momenta occurs. Such features, obtained 
considering (e~,2e _ ) electron-impact ionization within 
a quantum-mechanical S'-Matrix framework, are inter- 
preted in terms of the so-called quantum orbits, which 
can be directly associated to the trajectories of classical 
electrons. Both the asymmetry and the critical phase re- 
sult from the interplay between phase-space effects, and 
the probability that the first electron leaves its parent ion 
through tunneling ionization. The former and the latter, 
respectively, determine whether electron-impact ioniza- 
tion is classically allowed or forbidden, or the relevance 
of a set of orbits to the yield. 

The huge effects observed, namely the yields vanishing 
or appearing over extense and well-separated regions in 
the (pi||,2?2||) plane, are a due to a particular character- 
istic of the rescattering process in question, for which, in 
addition to a maximal, there is also a minimal classically 
allowed energy. In other words, electron- impact ioniza- 
tion is allowed only within confined phase-space regions, 
defined by the radius in Eq. (|ll)|l . By varying the driving- 
field intensity and the absolute phase adequately, this ra- 
dius can be forced to vanish, so that a whole phase-space 
region would become classically forbidden. Furthermore, 
a particular region can be made irrelevant due to a small 
tunneling probability for the first electron. This is a 
major advantage over other phenomena occurring in the 
context of strong-laser field matter interaction, such as 
above-threshold ionization and high-order harmonic gen- 
eration. For both phenomena, there are only maximal 
classically allowed energies, so that these effects do not 
occur. 

Apart from providing support for previous classical 
computations 29], the present results allow one to es- 
tablish a one-to-one correspondence between the classical 
and the quantum-mechanical approaches. 

Tunneling, for instance, is incorporated in the classical 
model using the quasi-static rate i)23|) which weights the 
electron counts. In the S'-Matrix computation, this pro- 
cess is directly related to the imaginary parts of the start 
times t'. Both the quasi-static rate and Im[t'] are some- 
how a measure of the relevance of a set of orbits. Indeed, 
dominant contributions always come from pairs of orbits 
for which Ea. (|23[l is larger, or Im[i'] are smaller than 
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those from the remaining pairs. Depending on the abso- 
lute phase, this may occur for a single set of orbits, whose 
contributions lie either in the negative or in the positive 
momentum regions, or there may be sets of orbits whose 
contributions compete. A critical phase characterizes a 
change in the dominant set of orbits. 

Furthermore, the phase-space effects which occur if 
electron-impact ionization becomes classically forbidden, 
i.e., if the radius in Eq. i|10|l collapses, are present in 
both frameworks in very similar, though not entirely 
identical, ways. In the classical computation, this would 
lead to vanishing yields, since the condition in the ar- 
gument of the S function in Eq. (|24|l would never be 
fulfilled. In the quantum-mechanical model, there would 
be exponentially decaying contributions throughout. If 
this radius does not collapse, the process will be allowed 
within a confined phase-space region. In the classical 
case, the start and return times coalesce at the bound- 
aries of this region, whereas, quantum-mechanically, this 
does not completely happen. 

The above-stated effects explain the minor discrep- 
ancies between Figs. 1 and 9. In particular, two 
types of discrepancies have been observed. First, for 
high driving-field intensities, the distributions obtained 
with the quantum-mechanical computation are slightly 
broader than those from the classical model. Second, 
in the quantum-mechanical framework, the distributions 
start to shift at a slightly smaller phase, as compared to 



the outcome of the classical simulation. The former dif- 
ferences were also present in NSDI with monochromatic 
driving fields, whereas the latter discrepancies are specific 
to asymmetric momentum distributions, which occur, for 
instance, in the context of few-cycle driving pulses. 

The agreement between both the quantum-mechanical 
and the classical computations go beyond the physical 
explanations for the asymmetry and the critical phase. 
Indeed, the momentum distributions computed with the 
one or the other method, as well as the predicted criti- 
cal phase and the interval in which the momenta change 
sign, are very similar, apart from minor differences near 
the classical boundaries. This is a concrete evidence that 
the effects reported in this paper are not rooted on a par- 
ticular model or framework, being, on the contrary, of a 
deeper physical nature. In fact, recently, a similar asym- 
metry has been observed in ongoing NSDI experiments 
with few-cycle driving pulses |4l| . 
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